Exponential Operators, Dobinski Relations and 
Summability 



P. Blasiak 1,2 , A. Gawron 1 , A. Horzela 1 , K. A. Penson 2 
and A. I. Solomon 2 3 

1 H.Niewodniczanski Institute of Nuclear Physics, Polish Academy of Sciences, 

ul. Eliasza-Radzikowskiego 152, PL 31342 Krakow, Poland 

2 Laboratoire de Physique Theorique de la Matiere Condensee, Universite 

P. et M. Curie, Tour 24 - 2e et., 4 Pl.Jussieu, F 75252 Paris Cedex 05, France 

3 Physics and Astronomy Department, The Open University, 

Milton Keynes MK7 6AA, United Kingdom 

E-mail: pawel . blasiakOif j . edu.pl, agnieszka.gawronOif j .edu.pl, 
andrzej .horzelaOif j .edu.pl, penson@lptl.jussieu.fr, 
a . i . solomonOopen .ac.uk 

Abstract. We investigate properties of exponential operators preserving the particle 
number using combinatorial methods developed in order to solve the boson normal 
ordering problem. In particular, we apply generalized Dobinski relations and methods 
of multivariate Bell polynomials which enable us to understand the meaning of 
perturbation-like expansions of exponential operators. Such expansions, obtained as 
formal power series, are everywhere divergent but the Pade summation method is 
shown to give results which very well agree with exact solutions got for simplified 
quantum models of the one mode bosonic systems. 
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Consider exponential operators of the form exp (— XTi, a (n)), where T~i a (n) = 
J2iLi a i n% is a polynomial of the number operator n (assumed to be bounded from 
below), Oj's play roles of coupling constants and A is an overall positive parameter. 
Calculation of the number state representation of exp (— XH a {nj) is straightforward 

N 

(l\exp(-\H a (h))\m) = exp (-A £ c^m 4 )^™, (1) 

i=\ 

as well as its (standard) coherent state representation (\z) = exp (— ■^j-) Y^=o l n )> 
h\n) = n\n), (n\n') = 5 nn ', a\z) = z\z) and (z'\z) = exp (— ~ (\z'\ 2 + \z\ 2 — 2z'z) s j) 

(z'\ex P (-XH a (n))\z) = (^} e -*'* £ L^L exp -A]>>F . (2) 

fc=0 fc! V i=l / 

For the simplest example of H a (n) = n the r.h.s of Eqn.Q is given in terms of the 

elementary function 

(z'|exp(-An)|z) = (z'\z)e~ z ' z f] ^p^exp(-A&;) 



a *» r; 7 (3) 

exp (z'z (e~ A - lJJ 



in which one recognizes the exponential generating function of the (exponential) Bell 
polynomials pQ 

oo \ k 

exp(o:(e A -i)) =]T-£(fc,x). (4) 
fc=o K - 

The Bell polynomials are well known from their applications in combinatorics [2|. They 
are defined as 

B(n,x) = Y,S(n,k)x k , (5) 

k=\ 

where S(n, k) denote the Stirling numbers of the second kind (positive integers which in 
enumerative combinatorics count the number of ways of putting n different objects into 
k identical containers leaving none container empty) whose analytic representation is 



^*) = ^E(*)(-i)*- i r 



r- (6) 

Particular values of the Bell polynomials B(n) = B(n, 1) are known as the Bell numbers 
and in enumerative combinatorics count the number of ways of putting n different 
objects into n identical containers some of which may be left empty This means that 
the Bell numbers give us the number of partitions of an n-element set. Expanding the 
l.h.s. of the Eqn.(jnj) as a power series in A and using the definition of the coherent states 
we arrive at 

oo \fc oo / ' \m oo \k 

(z'\z)e- Z ' Z ElfE [ -^T^ = (z'\z) £ ^B(k, z'z) (7) 

fc=0 m=0 " L - k=0 rh - 

from which one reads out the formula 

oo (J y \m 
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giving for z' 



1 the Dobirlski relation 

-.k 



Bin) 



-i 



m=0 



. m 
^ ml 



(9) 



Eqns. (JSJ) and Q connect sequences of polynomials B(k,z'z) or positive integers B(n) 
with sums of nontrivial series of fractions and allow to represent the Bell polynomials as 
the Stieltjes moments of an infinite sum of weighted 5-functions, called the Dime comb 



[■OC °° x m 

B(n,x) = e~ x / dyy 11 ^ — -5(y 
J o ml 



m) 



(10) 



Following the above considerations we can generalize our results to arbitrary 
exp (—XH a (n)) in which TL a {n) polynomially depend on h. For such a case we get 

:-A) fc 



z'\ exp {-XH a {n))\z) = (z'\z) £ B a (k, 



z'z) 



k=0 



kl 



(11) 



where generalized Bell polynomials B Q (m,z'z) are defined through generalized Stirling 
numbers of the second kind [3], jlj, 



mN 



S a (n, k) 




E S a (n, k) x k 



.l)k-j 



N 



E ai j- 



i=i 



(12) 



with j- — j(j — 1) . . . (J — I + 1) denoting the falling factorial and 6>i = Em=i m ) a m 
being the inverse Stirling transform of cti given by the Stirling numbers of the first kind 
s(n, k), Efc= m s ( n , k)S(k, m) = 5 nm . Generalized Dobinski relations read now j3], |3] 



B n 



m=0 



N 



m- 



ml 



(13) 



and, analogously to Eqn. pOjl . provide us with representation of B a (n,x) as moments 



„ oo 

B a (n,x)=e- X dyy n J2 



X 



ml 




where the domain of integration A = [inf (EfcLi a k m- 



■'>(/y-E''A '"-)• (14) 
e, oo). Eqn.((T3J), if put into 

Eqn. pi|) and changed the summation order, leads to Eqn.Q. Note that using the 
generalized Dobinski formula we give analytical meaning to the formal series (jlljl . As a 
rule these series are divergent because the coefficients B a (n, x) grow with n much faster 
than nl. Such an asymptotic behavior is seen from Eqns. (jl2j) - the latter imply that the 
numbers S a (n,k) include the standard Stirling numbers of the second kind S(nN,k) 
and, as a consequence, the polynomials B a (n,x) include polynomials B(nN,x). The 
n — > oo asymptotics of the standard Bell numbers Bin, 1) is B n ~ nl — cx p( cx pH.")) 1 ) 

J 1 V ' ' [r(n)]«+ 1 v /27rexp (r(n)) 

where r(n) ~ logn — log(logn) and it causes that the series (JTTj) are divergent for N > 2. 

For the toy model TC a (n) = n we were able to find the closed form of 
(z'\ exp (—An) 1 2) given in terms of elementary functions - i.e. we solved explicitly the 
normal ordering problem for such an operator jHj. If TC a (n) becomes a more complicated 
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polynomial then the problem complicates but it remains manageable and gives some 
insight into perturbation methods widely used in quantum mechanics and quantum field 
theory. Because in the following we are going to concentrate ourselves on the problems 
related to the coupling constant perturbation calculus treated with combinatorics-based 
methods we do not use the scheme leading to the generalized Bell polynomials but we 
will investigate the problem using methods of multivariate Bell polynomials, [2|, [Zj, 
still emphasizing the importance of the Dobihski-type relations. The multivariate Bell 
polynomials enable us to construct the Taylor-Maclaurin expansion of a composite 
function f(g{x)). To this end let us recall that for any f(x) = Z)^=i fnX n /n\ and 
g{ x ) — Hn=i 9nX n /n\ given as formal power series one gets 

oo I n \ n 

f(g( X )) = [f°g]( X ) = \J2 B nk(9l,92,---,9n-k+l)fk) — , (is) 

71=1 \fc = l / H - 

where the coefficients B nk are certain polynomials in the Taylor coefficients gi - namely 
the multivariate Bell polynomials - given by 

// ?tJ 

B nik ( 9l , . . . , g n ^ k+1 ) = J2 ~ : 9x92 ■ ■ ■ 9n-kTi , ( 16 ) 

n W (iO^'l 

where the summation " is over & h possible non- negative {z/j} being partitions of 

an integer n into sum of k integers, i.e. over {i^} being solutions to the equations 
YTj=\3 v j = n and J2]=i^j = k%. Eqns. (fT5|) and (fTBj) imply that the multivariate Bell 
polynomials satisfy, for a and b arbitrary constants, the homogeneity relation 

Bn^atfgx, ab 2 g 2 , ab n ~ k+1 g n -k+i) = a k b n B n>k (g l7 g 2 ..., g n -k+i) ■ (17) 

A particular case of the multivariate Bell polynomials are polynomials being coefficients 
of Taylor-Maclaurin expansions of exp(J2iLo ctiX 1 ). They generalize the standard Hermite 
polynomials and for some special cases have simple analytic forms [5] , JU| - among them 
the two variable Hermite-Kampe de Feriet polynomials H^ M \x,y): 

x n 



exp (gix + g M x M ) = HjV)(g 1} g M )— , 

[n/M] n-Mr-r (1 § ) 

HW( gi ,g M )=n\j: ' ;i <hl 



u (n — Mr)\r\ ' 
and the three- variable Hermite polynomials H n (a,b, c): 

exp (ax + bx 2 + cx 3 ) = H n (a, 6, c) — 



H n (a,b,c) =n\J2 



x' 
n\ 

c r Hl_ 3r (a,b) ^ 
, 1 1 K n - 3r)!r! 



In/3] Zu2 t„ u\ ( 19 ) 



\ This condition shows that the multivariate Bell polynomials are closely related to combinatorial 
numbers 
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As an illustration of the presented approach let us consider diagonal coherent state 
matrix element (z\ exp (— (g(£n) + G(C,n) 2 ))\z) . Expanding the exponential as power 
series in £, next using (j!7j) . Q18j) and definition of the coherent states we arrive at 

oo en 

(z\ exp (- ( 9 m + Gm 2 ))\z) = 1 + E itfH-g, ~G)B{n, k| 2 )^.(20) 

n=l U - 

The operational relation (x,y) = exp (y^j)x n enables us to rewrite the Eqn. (f2Tty 

as 

( z \e-('W+ G M a )\z) = exp(-G^)exp(|^| 2 (e-^-l)) 

= g^E5(/,N 2 )fr^(-^ (21) 
k=o "" i=o dp 

oo t2fc oo (—n£Y 

= E(-G) k ^ r EB(2k + i,\z\ 2 )^W-. 

Replacing B(2k + i, \z\ 2 ) by their moment representation (fTUj) and changing 
the integration and summation order we obtain the perturbation expansion of 
(z\ exp (— (g(£,n) + G(^n) 2 ^ \z) in terms of the power series in the coupling constant G 

(z\e-^ +G ^ 2 )\z) = e ^ 9 ^) g (- G )k(f) B (2k, |z| 2 e -^) (22) 

k=o K - 

Because of the asymptotic behavior B(2k, \z\ 2 e~ 9 ^) ~ (2k)\,k — > oo, this series, as 
well as the series (|20|) . both have zero radii of convergence being however asymptotic 
expansions of 

<z| exp (- U^) + G(en) 2 ))|.) = e-l^ 2 £ lfL e -(*(W^(W a ). (23) 

fc=o fc - 

In order to give them analytical meaning one has to use methods of generalized 
summation. Numerical check using the Pade method (see Fig. 1-2, below) shows that 
even low order approximants of series (|2Uj) and (J22j) give very good agreement with exact 
result calculated from ()23|) for G belonging to the domain much larger than the domain 
in which partial sums of both series give acceptable results. Moreover, comparing these 
results we in fact compare results given by the Pade method with exact solutions for 
simplified, nevertheless essentially quantum, model. This confirms practical utility of 
the Pade summation method applied to various perturbation expansions occurring in 
quantum physics even if we are unable to prove its applicability in a mathematically 
satisfactory way. It also confirms that generating functions obtained as solutions to 
the boson normal ordering problem and being in general divergent formal series may 
be interpreted as asymptotic expansions, resumed in such a generalized sense and, as a 
consequence, used in physical applications. 
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Figure 1. The sum of the series (12311 . g — £ = 1, (left plot) and its comparison with 
Pade approximants [3, 4] for the series (|20|l (dashed curve, worse approximation) and 
(|22|l (continuous curve, nondistinguishable from the plot of (|23|l ) for G = 1. 




Figure 2. Subdiagonal Pade approximants [3,4], g = £ = 1, for the series (|20|l (left 
plot) and (|22|l (right plot). It is seen that both approximations mimic the exact result 
(|23ll very well. It may be also checked that the results weakly depend on the order of 
approximation. 
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